Prominence seismology using the period ratio of transverse thread oscillations

The ratio of the period of the fundamental mode to that of the first overtone of kink oscillations, from here on the "period ratio", is a seismology tool that can be used to infer information about the spatial variation of density along solar magnetic flux tubes. The period ratio is 2 in longitudinally homogeneous thin tubes, but it differs from 2 due to longitudinal inhomogeneity. In this paper we investigate the period ratio in longitudinally inhomogeneous prominence threads and explore its implications for prominence seismology. We numerically solve the two-dimensional eigenvalue problem of kink oscillations in a model of a prominence thread. We take into account three nonuniform density profiles along the thread. In agreement with previous works that used simple piecewise constant density profiles, we find that the period ratio is larger than 2 in prominence threads. When the ratio of the central density to that at the footpoints is fixed, the period ratio depends strongly on the form of the density profile along the thread. The more concentrated the dense prominence plasma near the center of the tube, the larger the period ratio. However, the period ratio is found to be independent of the specific density profile when the spatially averaged density in the thread is the same for all the profiles. An empirical fit of the dependence of the period ratio on the average density is given and its use for prominence seismology is discussed.Comment: Accepted for publication in A&

Time damping of non-adiabatic magnetohydrodynamic waves in a partially ionized prominence plasma: Effect of helium

Prominences are partially ionized, magnetized plasmas embedded in the solar corona. Damped oscillations and propagating waves are commonly observed. These oscillations have been interpreted in terms of magnetohydrodynamic (MHD) waves. Ion-neutral collisions and non-adiabatic effects (radiation losses and thermal conduction) have been proposed as damping mechanisms. We study the effect of the presence of helium on the time damping of non-adiabatic MHD waves in a plasma composed by electrons, protons, neutral hydrogen, neutral helium (He I), and singly ionized helium (He II) in the single-fluid approximation. The dispersion relation of linear non-adiabatic MHD waves in a homogeneous, unbounded, and partially ionized prominence medium is derived. The period and the damping time of Alfven, slow, fast, and thermal waves are computed. A parametric study of the ratio of the damping time to the period with respect to the helium abundance is performed. The efficiency of ion-neutral collisions as well as thermal conduction is increased by the presence of helium. However, if realistic abundances of helium in prominences (~10%) are considered, this effect has a minor influence on the wave damping. The presence of helium can be safely neglected in studies of MHD waves in partially ionized prominence plasmas.Comment: Research note submitted in A&

The Thermal Instability of Solar Prominence Threads

The fine structure of solar prominences and filaments appears as thin and long threads in high-resolution images. In H-alpha observations of filaments, some threads can be observed for only 5 - 20 minutes before they seem to fade and eventually disappear, suggesting that these threads may have very short lifetimes. The presence of an instability might be the cause of this quick disappearance. Here, we study the thermal instability of prominence threads as an explanation of their sudden disappearance from H-alpha observations. We model a prominence thread as a magnetic tube with prominence conditions embedded in a coronal environment. We assume a variation of the physical properties in the transverse direction, so that the temperature and density continuously change from internal to external values in an inhomogeneous transitional layer representing the particular prominence-corona transition region (PCTR) of the thread. We use the nonadiabatic and resistive magnetohydrodynamic equations, which include terms due to thermal conduction parallel and perpendicular to the magnetic field, radiative losses, heating, and magnetic diffusion. We combine both analytical and numerical methods to study linear perturbations from the equilibrium state, focusing on unstable thermal solutions. We find that thermal modes are unstable in the PCTR for temperatures higher than 80,000 K, approximately. These modes are related to temperature disturbances that can lead to changes in the equilibrium due to rapid plasma heating or cooling. For typical prominence parameters, the instability time scale is of the order of a few minutes and is independent of the form of the temperature profile within the PCTR of the thread. This result indicates that thermal instability may play an important role for the short lifetimes of threads in the observations.Comment: Accepted for publication in Ap

Efficient implementation of a van der Waals density functional: Application to double-wall carbon nanotubes

We present an efficient implementation of the van der Waals density functional of Dion et al [Phys. Rev. Lett. 92, 246401 (2004)], which expresses the nonlocal correlation energy as a double spacial integral. We factorize the integration kernel and use fast Fourier transforms to evaluate the selfconsistent potential, total energy, and atomic forces, in N log(N) operations. The resulting overhead in total computational cost, over semilocal functionals, is very moderate for medium and large systems. We apply the method to calculate the binding energies and the barriers for relative translation and rotation in double-wall carbon nanotubes.Comment: 4 pages, 1 figure, 1 tabl

Kelvin-Helmholtz instability in partially ionized compressible plasmas

The Kelvin-Helmholtz Instability (KHI) has been observed in the solar atmosphere. Ion-neutral collisions may play a relevant role for the growth rate and evolution of the KHI in solar partially ionized plasmas as in, e.g., solar prominences. Here, we investigate the linear phase of the KHI at an interface between two partially ionized magnetized plasmas in the presence of a shear flow. The effects of ion-neutral collisions and compressibility are included in the analysis. We obtain the dispersion relation of the linear modes and perform parametric studies of the unstable solutions. We find that in the incompressible case the KHI is present for any velocity shear regardless the value of the collision frequency. In the compressible case, the domain of instability depends strongly on the plasma parameters, specially the collision frequency and the density contrast. For high collision frequencies and low density contrasts the KHI is present for super-Alfvenic velocity shear only. For high density contrasts the threshold velocity shear can be reduced to sub-Alfvenic values. For the particular case of turbulent plumes in prominences, we conclude that sub-Alfvenic flow velocities can trigger the KHI thanks to the ion-neutral coupling.Comment: Accepted for publication in Ap

Torsional Alfv\'en waves in solar partially ionized plasma: effects of neutral helium and stratification

Ion-neutral collisions may lead to the damping of Alfven waves in chromospheric and prominence plasmas. Neutral helium atoms enhance the damping in certain temperature interval, where the ratio of neutral helium and neutral hydrogen atoms is increased. Therefore, the height-dependence of ionization degrees of hydrogen and helium may influence the damping rate of Alfven waves. We aim to study the effect of neutral helium in the damping of Alfven waves in stratified partially ionized plasma of the solar chromosphere. We consider a magnetic flux tube, which is expanded up to 1000 km height and then becomes vertical due to merging with neighboring tubes, and study the dynamics of linear torsional Alfven waves in the presence of neutral hydrogen and neutral helium atoms. We start with three-fluid description of plasma and consequently derive single-fluid magnetohydrodynamic (MHD) equations for torsional Alfven waves. Thin flux tube approximation allows to obtain the dispersion relation of the waves in the lower part of tubes, while the spatial dependence of steady-state Alfven waves is governed by Bessel type equation in the upper part of tubes. Consecutive derivation of single-fluid MHD equations results in a new Cowling diffusion coefficient in the presence of neutral helium which is different from previously used one. We found that shorter-period (< 5 s) torsional Alfven waves damp quickly in the chromospheric network due to ion-neutral collision. On the other hand, longer-period (> 5 s) waves do not reach the transition region as they become evanescent at lower heights in the network cores. Propagation of torsional Alfven waves through the chromosphere into the solar corona should be considered with caution: low-frequency waves are evanescent due to the stratification, while high-frequency waves are damped due to ion neutral collisions.Comment: 9 pages, 7 figures (accepted in A&A